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Abstract. We study the incompressible limit of the compressible non- 
iscntropic ideal magnctohydrodynamic equations with general initial data in 
the whole space R d (d = 2,3). We first establish the existence of classic solu- 
tions on a time interval independent of the Mach number. Then, by deriving 
uniform a priori estimates, we obtain the convergence of the solution to that 
of the incompressible magnctohydrodynamic equations as the Mach number 
tends to zero. 



1. Introduction 

This paper is concerned with the incompressible limit to the compressible ideal 
non-isentropic magnctohydrodynamic (MHD) equations with general initial data in 
the whole space M. d (d = 2, 3), which take the following form 

dtp + div(pu) = 0, (1.1) 

dt(pu) + div (pu <g> u) + Vp = (V x H) x H, (1.2) 

d,H- V x (u x H) = 0, divH = 0, (1.3) 

d t £ + div (u(f ' + p)) = div((u x H) x H). (1.4) 

Here p denotes the density, u € R d the velocity, and H G K d the magnetic field, 
respectively; £ is the total energy given by £ = £' + |H| 2 /2 and £' = p (e + |u| 2 /2) 
with e being the internal energy, p|u| 2 /2 the kinetic energy, and |H| 2 /2 the magnetic 
energy. The equations of state p = p(p, 9) and e = e(p, 9) relate the pressure p and 
the internal energy e to the density p and the temperature 9 of the flow. 

For smooth solutions to the system (|1.1|) - (|1.4|) . we can rewrite the total energy 
equation (|1.4[) in the form of the internal energy. In fact, multiplying (|1.2[) by u and 
(|1.3[) by H respectively, and summing the resulting equations together, we obtain 



i|H| 2 ) + ldiv(p|u| 2 u) + Vp • u 



= (V x H) x H -u + V x (u x H) • H. (1.5) 

Using the identities 

div(H x (V x H)) = |V x H| 2 - V x (V x H) • H, (1.6) 

div((u x H) x H) = (V x H) x H • u + V x (u x H) • H, (1.7) 
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and subtracting (|1.5[) from (|1.4p . wc thus obtain the internal energy equation 

dt(pe) + div(pue) + (divu)p = 0. (1.8) 
Using the Gibbs relation 

6dS = de + pd(^J , (1.9) 

we can further replace the equation (|1.8[) by 

d t {pS) + div(pSu) = 0, (1.10) 

where S denotes the entropy. 

Now, as in |26j . we reconsider the equations of state as functions of S and p, 
i.e., p = R(S,p) for some positive smooth function R defined for all S and p > 
satisfying dR/dp > 0. For instance, we have p = p 1 ' 7 e -fif ' r w jth 7 > 1 for a 
polytropic gas. Then, by utilizing (|1.1[) together with the constraint divH = 0, the 
system (|l.ip - (|1.3p and (| 1 . 10[) can be rewritten as 

A(S,p)(d t p + (u • V)p) + divu = 0, (1.11) 
R(S,p)(d t u + (u • V)u) + Vp = (V x H) x H, (1.12) 
d t H - curl (u x H) = 0, divH = 0, (1.13) 
d t S + (u- V)S = 0, (1.14) 

where 

Considering the physical explanation of the incompressible limit, we introduce 
the dimensionless parameter e, the Mach number, and make the following changes 
of variables: 

p{x,t) =p t (x,et) =pe^^ x ^\ S(x,t) = S* £ (a;,ei), 
u(x, t) = eu e (x, et), H(x, t) = eH e (x, et), 

for some positive constant p. Under these changes of variables, the system (jl.lll) 
(|1.14|) becomes 

a(S e , eq e ){dtq e + (u e • V)g 6 ) + -Uivu £ = 0, (1.15) 

r(S e ,eq e )(d t u e + (u e ■ V)u e ) + -Vq e = (curlH 6 ) x H £ , (1.16) 

e 

d t U e - curl(u e x H e ) = 0, divH e = 0, (1.17) 

d t S e + (u e • V)5 £ = 0, (1.18) 



where we have used the abbreviations 

a(S e ,eq e ) := A(S e ,pe e,1 ')pe eqE 

R(S^pe^) 
r(S,eq):= ^ ■ 



pe eq dR(S e ,s) 



R(S e ,pe e i € ) ds 
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Formally, we obtain from (|1.15|) and ([1.16)1 that Vg e — > and divu £ — > as 
e — > 0. Applying the operator curl to ([1.16)1 . using the fact that curlV = 0, and 
letting e — > 0, we obtain that 

curl (r(S, Q)(d t v + v • Vv) - (curlH) x H) = 0, 

where we have assumed that (S e , q e , u £ , H £ ) and r(S e ,eq e ) converge to (5, 0, v, H) 
and r(S,0) in some sense, respectively. Finally, applying the identity 

curl (u x H) = u(divH) - H(divu) + (H • V)u - (u • V)H, (1.19) 

we expect to get, as e — > 0, the following incompressible ideal non-iscntropic MHD 
equations 

r(S,0)(d t v + (v • V)v) - (curlH) x H + Vtt = 0, (1.20) 
3 t H + (v • V)H - (H • V)v = 0, (1.21) 
d t S + (v • V)S = 0, (1.22) 
divv = 0, divH = (1.23) 

for some function it. 

The aim of this paper is to establish the above limit process rigorously in the 
whole space M d . 

Before stating our main result, we briefly review the previous related works. 
We begin with the results for the Euler and Navier-Stokes equations. For well- 
prepared initial data, Schochet [32] obtained the convergence of the compressible 
non-isentropic Euler equations to the incompressible non-isentropic Euler equations 
in a bounded domain for local smooth solutions. For general initial data, Metivier 
and Schochet [26] proved rigorously the incompressible limit of the compressible 
non-isentropic Euler equations in the whole space K d . There are two key points in 
the article First, they obtained the uniform estimates in Sobolev norms for 

the acoustic component of the solutions, which are propagated by a wave equation 
with unknown variable coefficients. Second, they proved that the local energy of the 
acoustic wave decays to zero in the whole space case. This approach was extended 
to the non-isentropic Euler equations in an exterior domain and the full Navier- 
Stokes equations in the whole space by Alazard in [JJ and [2J, respectively, and to 
the dispersive Navier-Stokes equations by Levermore, Sun and Trivisa )24j . For the 
spatially periodic case, Metivier and Schochet [27] showed the incompressible limit 
of the one-dimensional non-iscntropic Euler equations with general data. Compared 
to the non-isentropic case, the treatment of the propagation of oscillations in the 
isentropic case is simpler and there is an extensive literature on this topic. For 
example, see Ukai [34], Asano [3], Desjardins and Grenier [7] in the whole space; 
Isozaki [14l[T5] in an exterior domain; Iguchi [13] in the half space; Schochet [31] 
and Gallagher [Til] in a periodic domain; and Lions and Masmoudi [2S], and Des- 
jardins, ct al. in a bounded domain. Recently, Jiang and Ou [SO] investigated 
the incompressible limit of the non-iscntropic Navier-Stokes equations with zero 
heat conductivity and well-prepared initial data in three-dimensional bounded do- 
mains. The justification of the incompressible limit of the non-isentropic Euler or 
Navier-Stokes equations with general initial data in a bounded domain or a multi- 
dimensional periodic domain is still open. The interested reader can refer to [5] on 
formal computations for the case of viscous polytropic gases and [U[27] on some 
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analysis for the non-isentropic Euler equations in a multi-dimensional periodic do- 
main. For more results on the incompressible limit of the Euler and Navier-Stokes 
equations, please see the monograph [S] and the survey articles pl l5Ull55] . 

For the isentropic compressible MHD equations, the justification of the low Mach 
limit has been given in several aspects. In [21] . Klainerman and Majda first studied 
the incompressible limit of the isentropic compressible ideal MHD equations in the 
spatially periodic case with well-prepared initial data. Recently, the incompressible 
limit of the isentropic viscous (including both viscosity and magnetic diffusivity) 
of compressible MHD equations with general data was studied in [T^IIISIIIZ] • In 
[12] , Hu and Wang obtained the convergence of weak solutions of the compressible 
viscous MHD equations in bounded, spatially periodic domains and the whole space, 
respectively. In [16] . the authors employed the modulated energy method to verify 
the limit of weak solutions of the compressible MHD equations in the torus to 
the strong solution of the incompressible viscous or partial viscous MHD equations 
(the shear viscosity coefficient is zero but the magnetic diffusion coefficient is a 
positive constant). In |17j . the authors obtained the convergence of weak solutions 
of the viscous compressible MHD equations to the strong solution of the ideal 
incompressible MHD equations in the whole space by using the dispersion property 
of the wave equation if both shear viscosity and magnetic diffusion coefficients go 
to zero. 

For the full compressible MHD equations, the incompressible limit in the frame- 
work of the so-called variational solutions was studied in [22j[23l[28]. Recently, 
the authors [18] justified rigourously the low Mach number limit of classical solu- 
tions to the ideal or full compressible non-isentropic MHD equations with small 
entropy or temperature variations. When the heat conductivity and large temper- 
ature variations are present, the low Mach number limit for the full compressible 
non-isentropic MHD equations was shown in [19| . We emphasize here that the 
arguments in |19] are completely different from the present paper (at least in the 
derivation of the uniform estimates), and depend essentially on the positivity of 
fluid viscosities, magnetic diffusivity, and heat conductivity coefficients. 

As aforementioned, in this paper we want to establish rigorously the limit as 
e — > for the system (|1.15[) - (|1.18[) . In this case, there arc additional magnetic field 
terms appearing in the momentum equations, compared with non-isentropic Euler 
equations. When applying the non-standard derivative operators to the momentum 
equations, we have to cancel the the most-differentiated terms to close our estimates 
with the help of the magnetic field equations. However, it is difficult to show this 
cancelation due to the strong coupling of hydrodynamic motion and magnetic field. 
To surround such difficulties, new ideas and techniques are required in the present 
paper, and we shall explain these briefly after stating our result. 

Now, we supplement the system (|1.15[) - (|1.18j) with initial conditions 

(SW,u e ,H e )| t=0 = (S e ,q e ,ul,Hl). (1.24) 
Our main result thus reads as follows. 

Theorem 1.1. Suppose that the initial data .Sq , gQ , Uq , Hq satisfy 

|[(S§,g§ > u§,H§)|| H « (R , ) <M (1.25) 

for some constant Mo > 0. Then there exist constants T > and e$ G (0, 1) 
such that for any e € (0, eg], the Cauchy problem (|1.15j) - (|1.18|) . (]1.24j) has a unique 
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solution (S e , q e , u e , H c ) G C*°([0, T], H 4 (BL d )), and there exists a positive constant 
N , depending only on T , eo, and Mo, such that 

\\(S\q\u\W)(t)\\ Hi{Kd) <N, Vfe[0,T]. (1.26) 

Furthermore, suppose that the initial data (Sq,Uq,H ) converge in H 4 (M. d ) to (So, 
w , H ) as e — > 0, and that there exist positive constants S_, Nq and 8, such that Sg 
satisfies 

\Sq(x)-S\ < Nolx^ 1 - 8 , \VS e (x)\ <N \x\- 2 - s , (1.27) 
then the sequence of solutions (S e , q e , u £ , H e ) converges weakly in L°°(0, T; H 4: (M d )) 
and strongly in L 2 (Q,T; H^ oc (M. d )) for all s' < 4 to a limit (S, 0,v,H), where 
(#,v,H) is the unique solution in C([0, T], H 4 (R d )) of (fTT^Uj) - pT^H]) with initial 
data (S,v, H)| t=0 = (So,Wo,Ho), where wo € H 4 (M. d ) is determined by 

divw o = 0, curl (r{S , 0)w ) = curl (r(5 ,0)v ), r(S , 0) := lim r(S%, 0). (1.28) 
T/ie function tt £ C*([0,T] x M d ) sote^es Vtt e C([0, T], H 3 (R d )). 

We briefly describe the strategy of the proof. The proof of Theorem 11.11 in- 
cludes two main steps: uniform estimates of the solutions, and the convergence 
from the original scaling equations to the limiting ones. Once we have established 
the uniform estimates (|1.26|) in Theorem 11.11 the convergence of solutions is easily 
obtained by using the local energy decay theorem for fast waves in the whole space 
shown by Metivier and Schochet in |26j . Thus, the main task in the present paper 
is to obtain the uniform estimates (|1.26p . For this purpose, we shall modify the 
approach developed in |26| . In the process of deriving uniform estimates, when we 
perform the operator ({E}~ 1 L(d x ))' T (c = 1,2,3,4) to the continuity and momen- 
tum equations, or the operator curl to the momentum equations, one order higher 
spatial derivatives arise for the magnetic field, and the key point in the derivation 
of the uniform estimates is to cancel these troublesome magnetic terms in the fluid 
equations by the corresponding terms in the magnetic field equations. Only those 
magnetic terms of highest order derivatives are needed to be controlled (canceled) . 
To achieve this cancellation, we restrict the Sobolev index s to be even, so that the 
highest order derivatives applied to the momentum equations are not intertwined 
with the pressure equation, and we then apply the same highest order derivative 
operators to the magnetic field equations. We remark that the divergence-free 
property of the magnetic field plays an important role in the process of the uniform 
estimates. 

Remark 1.1. In the proof of Theorem ll.il it is crucial to appropriately choose the 
index s of the Sobolev space H s to be even (here we take s = 4 for simplicity). The 
result is still valid for any even integer s > d/2+1. In contrast, for the non-isentropic 
Euler equations, the same result holds in H s for any integer s > d/2 + 1 (see [26]). 
The restriction on the index is due to the strong coupling of hydrodynamic motion 
and magnetic field. 

Remark 1.2. The uniform estimates obtained in this paper are still valid when we 
consider our problem in a torus T d (d = 2, 3) . However, it is an open problem to 
prove rigorously the incompressible limit of the compressible non-isentropic ideal 
MHD equations in higher dimension T d due to the loss of dispersive estimates 
for the wave equations, see |27| for some discussions on the non-isentropic Euler 
equations. 
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Remark 1.3. We point out that our arguments in this paper can be modified slightly 
to the case of the non-isentropic MHD equations with nonnegative fluid viscosities 
or magnetic diffusivity but without heat conductivity. In this case we do not require 
that the index s of the Sobolcv space H s must be even. The reason is that we can use 
the diffusion term (terms) to control the highest troublesome terms. For example, 
considering the situation that the magnetic diffusivity is zero and the fluid viscosity 
coefficients are positive, we can define 

M e (T) :=K(T) 2 + Aul^dr 
J o 

with 

Ne{T):= sup ||(5 e , g e ,u e ,H e )(i)|| s , s > d/2 + 2, 

te[o,T] 

and employ the same arguments with slight modifications to obtain the same result 
as Theorem 12. II 

Remark 1.4. It seems that the proof of the Theorem ll.ll can not be extended directly 
to the full compressible MHD equations where the heat conductivity is positive. In 
this situation we can not write the total energy equation as a transport equation 
of entropy, which plays a key role in the proof. Instead, it is more convenient to 
use the temperature as an unknown in the total energy equation when the heat 
conductivity is positive, see [19] for more details. 

Notions. We denote by (■, ■} the standard inner product in L 2 (M. d ) with norm 
Il/H 2 = (/,/} and by H k the usual Sobolev space W k ' 2 (R d ) with norm || ■ \\ k . 
In particular, j| ■ ||o = || ■ ||. The notation . . . , Ak)\\ means the summation 

of || (i = 1, ■•• ,k), and it also applies to other norms. For a multi-index 
a = (ai, . . . , ad), we denote d a = c^ 1 . . . and |a| = \ct\ \ + ■ ■ ■ + \ctd\- We shall 
omit the spatial domain M d in integrals for convenience. We use the symbols K 
or Co to denote generic positive constants, and C(-) and C(-) to denote smooth 
functions of their arguments which may vary from line to line. 

This paper is arranged as follows. In Section 2, we establish the uniform bound- 
eness of the solutions and prove the existence part of Theorem 11.11 In Section 3, 
we use the decay of the local energy to the acoustic wave equations to prove the 
convergent part of Theorem 11.11 

2. Uniform estimates 

Throughout this section e > 0, will be fixed, and the solution (S e , q e , u e , H c ) will 
be denoted by (S, q, u, H) and the corresponding superscript e in other notations 
are omitted for simplicity of presentation. 

In view of [26] and the classical local existence result in [25] for hyperbolic 
systems, we see that the key point in the proof of the existence part of Theorem ll.il 
is to establish the uniform estimate (|1.26[) . which can be deduced from the following 
a priori estimate. 

Theorem 2.1. Let (S,q,u,H) £ C([0, T], iJ 4 (M d )) be a solution to ([TTT5]) - ([i~18]) . 
(|1.24|) . Then there exist constants Co > 0, < eo < 1 and an increasing function 
C(-) from [0, oo) to [0, oo), such that for all e £ (0, eo] and t £ [0, T], 

M(T) < Co + (T + e)C(M(T)), (2.1) 
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where 

M(T):= sup ||(S > g,u J H)(t)|| 4 . (2.2) 

te[o,T] 

The remainder of this section is devoted to establishing (|2.1|) . In the calculations 
that follow, we always suppose that the assumptions in Theorem 11.11 hold. We 
consider a solution (S,q,u,H) £ C([0, T], H 4 (R d )) to the problem ([TT5]) - (p38]) . 
(11.241) with the initial data satisfying (|1.25p . 

First, we have the following estimate of the entropy S, which was obtained in [26] . 

Lemma 2.2. There exist a constant Cq > and a function C(-), independent of 
e, such that for all e £ (0, 1] and t 6 [0, T], 

\\S(t)\\l<Co+tC(M(T)). (2.3) 

The following L 2 -bound of (q, u, H) can be obtained directly using the energy 
method due to the skew-symmetry of the singular term in the system and the 
special structure of coupling between the magnetic field and fluid velocity. This 
i 2 -bound is very important in our arguments, since it is the starting point to the 
following induction analysis used to get the desired Sobolcv estimates. 

Lemma 2.3. There exist a constant Cq > and an increasing function C(-) inde- 
pendent of e, such that for all e € (0, 1] and t £ [0, T\, 

||(g,u,H)(t)|| 2 <C + tC(M(T)). (2.4) 

Proof. Multiplying (|1.15[) by q, (|1.16l) by u, and (|1.17[) by H, respectively, integrat- 
ing over R d , and adding the resulting equations together, we obtain 

(ad t q,q) + (rd t u,u) + (d t U,H) 

+ (a(u • V)q, q) + (r(u • V)u, u) 

= J [(curl H) x H] • u dx + J curl (u x H) • H dx. (2.5) 

Here the singular terms involving 1/e are canceled. Using the identity (|1.7|) and 
integrating by parts, we immediately obtain 

J [(curlH) x H] ■ udx + J curl (uxH)-Hdi = 0. 

In view of the positivity and smoothness of a(S, eq) and r(S, eq), we get directly 
from (| 1 . 1 5|) . (|1.18|) and the well-known nonlinear estimates (e.g., see [IT]): 

||/(«)|U < C{\\u\\ a )\\u\\ a , u e H°(R d ), a > d/2, /(■) smooth with /(0) = 0, 

that 

\\d t S\\ 3 < C{M(T)), \\ed t qh < C{M(T)), (2.6) 

while by the Sobolev embedding theorem, we find that 

\\(d t a,d t r)\\ L - < \\(d t a,d t r)h < C(M(T)). (2.7) 

By the definition of M(T) and the Sobolev embedding theorem, it is easy to see 
that 

||(Va,Vr)|| L ~ <C(M(T)). 
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Thus, from (|2.5[) we get that 

(aq, q)+{m, u) + (H, H) < {(aq, q) + (ru, u) + (H, H)}| t=Q 

+ C(M(T)) f {|g(r)| 2 + |u(r)| 2 + |H £ (r)| 2 }dr. (2.8) 
Jo 

Moreover, we have 

||gir + ||u|| 2 <||(o)- 1 |U-(a ?J g> + ||(r)- 1 || i -<ru,u) 
<C ((ag,<z) + (ru,u}), 

since a and r are uniformly bounded away from zero. Applying Gronwall's Lemma 
to (|2.8p . we conclude that 

\\(q,u,H)(t)\\ 2 < Co\\(qo,n ,U Q )\\ 2 cxp{tC(M(T))}. 

Therefore, the estimate (|2.4j) follows from an elementary inequality 

e ct < 1 + Ct, 0<t< T , (2.9) 

where To is some fixed constant. □ 

To derive the desired higher order estimates, we shall adapt and modify the 
techniques developed in |26j . Set 

E (^)=("< s ^> r{s ° eq)ld ), k*>-(° r). »=(;)• 

where 1^ denotes the d x d unit matrix. 

Let L E {d x ) = {EiS^eq)}- 1 ^^ and r (5) = r(S,0). Note that r (5) is 
smooth, positive, and bounded away from zero with respect to each e. First, using 
Lemma 12.21 and employing the same analysis as in [26) , we have 

Lemma 2.4. There exist constants C\ > 0, K > ; and an increasing function 
C(-), depending only on Mq, such that for all a G {1,2,3,4}, e € (0,1] and t £ 

[o,n 

||U|| ff < K\\L{d x )V\\ a -i + CdlcurUrouJIU-i + ||U|| ff _i) (2.10) 

and 

||U|U < C{||{L B (a 2 0ru|| o + HcurUrou)!!^! + ||U|| ff _i}, (2.11) 
where C := d + tC(M(T)) + eC(M(T)). 

We remark that the inequalities (|2.10[) and (|2.1ip are similar to the well known 
Helmholtz decomposition, and the estimate on HS 1 ^)!! 2 in Lemma T2.2I plays a key 
role in the proof of Lemma 12.41 

Our next task is to bound ||{i£;(9 a; )} <T U||o and ||curl (roil) \\a-i by induction 
arguments. We first show the following estimate. 

Lemma 2.5. There exists a constant Co > and an increasing function C(-) from 
[0, oo ) to [0, oo ), independent of e, such that for all < a < 3, all e G (0,1] and 
t € [0, T], it holds that 

J E\{L E {d x )Y\5{t)\ 2 Ax < C +tC(M(T)), (2.12) 
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and for a = 4, 

J E\{L E {d x )YV(t)\ 2 dx 

<C +tC(M(T)) + 2 y"(ar)- 4 (V(A 2 H)H)VAdivudxdr. (2.13) 

Proof. Let U CT := {L E (d x )y\J, a E {0, . . . , 4}. For simplicity, we set M := M(T), 
and E := E(S,eq). The case k = is an immediate consequence of Lemma [2.31 
It is easy to verify that the operator LE(d x ) is bounded from H k to H k ~ 1 for 
k 6 {1, . . . , 5}. Note that the equations (|1.15[) . (|1.16[) can be written as 

(d t +u- \7)\J + -E- 1 L(d x )U = E-^-J (2.14) 
e 

with 



(curlH) x H 

For a > 1, we commute the operator {Lb} 17 with (|2.14l) and multiply the result- 
ing system by E to infer that 

E(d t + u • V)U CT + -L(d x )V a = E(f a + g CT ), (2.15) 
e 

where 

f CT := [d t + u- V, {L E y\V, g CT := {L J5 } CT ( J B- 1 J). 

Multiplying (|2.15p by U„ and integrating over (0, t) x R d with t <T, noticing 
that the singular terms cancel out since L(d x ) is skew-adjoint, we use the inequal- 
ities (|2.6|) and (|2.7|) . and Cauchy-Schwarz's inequality to deduce that 

±(E(t)U a (t),\J a (t)) <i( J B(0)U CT (0),U CT (0))+C(X)^||U ff MI| 2 dr 

+ /"||f<x(T)|| 2 dT+ /" / S(g CT U ff )(r)dr. (2.16) 
Jo Jo J 

For a < 3, the estimate (|2.12j) is obtained by taking a similar analysis to that of 
Lemma 2.4 in [26] . 

Next we mainly deal with the case a = 4. Following the proof process of Lemma 
2.4 in [21], we obtain that 

\m)\\<C(M(t)). (2.17) 

Now we estimate the nonlinear term in (|2.16[) involving g a . A straightforward 
computation implies that 



{LiL 2 } 2 Liu 
{L 2 Li} — L 2 q 



{L 2 L x Y l2 u 



if k is odd; 



if k is even, 



where 

Li := a _1 div, L 2 := r _1 V. 
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Therefore for k = 4, we can replace the operator {Ls(9 x )} 4 by 
{L E (d x )} 4 = 1 

(ar)- 2 A 2 







/ div 




(i 


i) 




v)] 



(ar)- 2 (Vdiv) 2 J 

where £ is a matrix whose elements are the summation of the lower order differ- 
ential operators up to order 3. Thus by virtue of Cauchy-Schwarz's and Sobolev's 
inequalities, and (|2.7p . we infer that 

J J Eg 4 XJ 4 dxdT = J £:(L s ) 4 (£:- 1 J)U4da;dT 

< J J (ar)- 4 [(Vdiv) 2 (curlH x H)](Vdiv) 2 uda;dr + C(M) 

=:h(t) + C(M). 
To control Ii(t), we recall the basic vector identities: 

div(a x b) = b • curia — a • curlb, curl curl a = V div a — Aa, 
and the fact that divH = to deduce 

h(t) < J J (ar)~ 4 (V(A 2 H)H)VAdivudxdr + C(7W). 



Thus, we have 

J J £;g 4 U 4 dxdT < J J (ar)- 4 (V(A 2 H)H)VAdivud.TdT + C(7W). (2.18) 

Finally, ([2.13|l follows from the above estimates (|2.16jl - (|2.18j) and the positivity 
of E. □ 

Next, we derive an estimate for ||curl(rou)|| ff _i. Define 

« s - e " : = 1 -^)- <" 9 » 

Hereafter we denote ro(t) := ro(S(t)) and f(t) := f(S(t),eq(t)) for notational 
simplicity. 

One can factor out eq in fit). In fact, using Taylor's expansion, one obtains that 
there exists a smooth function g(t), such that 

f(t) = eg(t) := eg(S(t), eq(t)), \\g(t)\\ 4 < C(M(T)). (2.20) 

Since 

d t S + (u ■ V)S = 0, 
the momentum equations (|1 . 16[) are equivalent to 

[d t + (u ■ V)](r u) + -Vq = gVq + (1 - e 5 )(curlH) x H. (2.21) 
e 

We perform the operator curl to the system (|2.2ip to obtain that 

[$ + (u-V)](curl(r u)) 
= [u • V,curl](r u) + [cml,g]Vq + curl[(l - eg) (curl H) x H], (2.22) 
and we can show 
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Lemma 2.6. There exist constants Co > ; < eo < 1 and an increasing function 
C(-) from [0,oo) to [0, oo), such that for all e £ (0, eo] and t € [0,T], it holds that 

|| {curl (r u) || \ < C +tC(M(T)), (2.23) 

and for \a\ = 3, 

2 (9 a curl (r u)(t), F<9 Q curl (r u)(t)> < C + tC(7W(T)) 

"1=3 

+ 2^] /" J{ar)- 4 (H- V)d a Ud a {Au)dxdT, (2.24) 

|a|=3 

where we denote by F := (ar)~ 4 (ro(l — eg)) -1 > 0. 

Proof. The estimate (|2.23|) directly follows from an energy estimate performed on 
(|2.22[) , and hence we omit the details here. Here we give a proof of (|2.24p only. 

Set M := M(T), and lu = curl(r u). By the defination of f(S,eq) and (|230| . 
there exists an eo € (0, 1], such that 1 — eg is bounded from below for e G (0, eo]. 
Taking d™ (\a\ = 3) to (|2.22|) . multiplying the resulting equations by Fd"u, and 
integrating over (0, t) x M. d with t < T, we obtain that 

J(S a w(t),F9 a o;(t)) <i(c) a w(0),Fa Q w(0)) +C(M) / ||a Q w(r)|| 2 dr 
2 2 Jo 

- / <(u e • V)d a u,Fd a uj)(T)dT 
Jo 

+ f ([u.V,t)>,F9 a W )(r)dr 
Jo 

+ [ (d a {[cur\,g e }\7q},Fd a Gj)(T)dT 
Jo 

+ [ (d a {[u e ■ V,curl](r u)}, J F i a a w)(r)dr 
Jo 

+ [ (<9 Q {curl[(l -e 3 )(curlH) x H]},F5 Q w)(r)dr 
Jo 

= : hd a u(0) 7 Fd a u>(0)} + f ][>(t). (2.25) 

We have to estimate the terms Ii(r) (2 < i < 6) on the right-hand side of (|2.25[) . 
Applying integrations by parts, we have 

h{r)<l F\d a ui\ 2 divudx + C(M) 

< ||divu(r)||io= \\d a uj(T)\\ 2 + C(M) < C(M), (2.26) 
while for the term ^(r), an application of Cauchy-Schwarz's inequality gives 
\h(r)\ < CWFd^W ||h a (r)||, h a (r) := [u • V,9>. 

The commutator h a is a sum of terms d^udjui with multi- indices /3 and 7 satisfying 
101 + h\ < 4, 1^1 > 0, and M > 0. Thus, 

||h Q (r)|| <C(X), 
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where the the following nonlinear Sobolev inequality has been used (see [H]): For 
all a = (ai, • • • ,ay), a > 0, and f,g 6 i/' £+ ' T (R' :i ), |a| = fe, it holds that 

||[/ j a a ] fl || ff <<7o(||/|| W ri.c.|| 5 || <r+fc _ 1 + ||/||« r+fc || fl || £ o.). 

Hence, we have 

\I 3 (r)\<C(M) + \\d a u J (r)r. (2.27) 

Noting that ([curl , g] a )ij = aid Xj g — ajd Xi g for a = (ai,--- ,a,d), using the 
estimate (|2.20[) . and the inequality ||mj|[ (T _fc_j < lf||u|| (T _fc||t>|| (r _j for k > 0, I > 0, 
k + 1 < a and a > d/2 (see [H]), we can control the term Ia(t) as follows. 

\h(r)\ < ||d Q {[curl,. ]Vg}|| WFd^W 

<K\\[cuT\,g]Vq\\ 3 \\Fd a Lj\\ 

<if||V fl (r)||3||Vg(r)||3||F5°a;|| 

<C(M) + \\d a uj(T)\\ 2 . (2.28) 

Similarly, the term Is(t) can be bounded as follows. 

|/ 6 (t)| < A-||d Q {[u-V,curl](r u)}|| \\F8 a oj\\ 

< Al[u.V,curl](r u e )|| 3 ||Fa Q W || 

<K\\[n j ,cml]d Xj {r u)\\ 3 \\Fd a u\\ 

<C(M) + \\d a Lu(T)\\ 2 . (2.29) 

For the term 1$ (r) , by virtue of the formula curl curl a = V div a — Aa, we 
integrate by parts to see that 

h(r) = ( ro ( ( a i r r e g) a " curl {( 1 - e 9)(( n ■ V ) H - V(ffi)}, 9«curl(r u)) 

< ~( 4 - eg)H ■ V)H},a a curlcurl(r u)\ + C(.M) 

^o(l-e.g) / 

< -((ar)- 4 9 a [(H ■ V)]H, <9 Q (Vdivu - Au)) + C(M) 

< ((ar)- 4 (H • V)<9 Q H, <9 a Au) + C(M), 

where we have used the fact that divH = 0. 

Inserting the estimates for Ij(r) (2 < j < 6) into (|2.25[) and summing over 
|a| = 3, we obtain (|2.24[) . The lemma is proved. □ 

We proceed to estimating the magnetic field. 

Lemma 2.7. There exist a constant Cq > and an increasing function C(-) from 
[0, oo) to [0, oo), such that for \a\ < 2, all e S (0, 1] and t £ [0, T], u;e /iat;e 

53 /(ar)- 4 |5 a curlH| 2 dx<Co + iC(A4(T)) (2.30) 

|a|<2 

and for \a\ = 3, 



5] /"(arr 4 |d a curlH| 2 d.T 
<C + tC{M{T)) - 2 53 j£ y M _4 (H • V)o> Q Ha Q (Au) dxdr 



|a|=3- 
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-2 J y"(ar)- 4 (VA 2 HH)VAdivuda:dT. (2.31) 

Proof. Set M := M.(T). We apply the operator curl to (|1.17j) and use the vector 
identity (|1.19|) to obtain that 

<9 t (curlH) + u- V(curlH) = -[curl,u] • VH + curl((H • V)u - Hdivu). (2.32) 

Taking <9 Q (|a| < 3) to (|2.32p . and multiplying the resulting equations by (ar)~ 4 x 
9 a curlH, and integrating over (0,i) x M. d with t < T, we find that 

((ar)- 4 d Q curl H, <9 Q curl H) 
= {({a(t)r(t)}- 4 a Q eurlH(0,a Q curlH(t))}| t=0 
+ (a t {(ar)" 4 }9 Q curlH,9 Q curlH) 

-2 J J {ar)- 4 d a [(u- V)curl H)<9 a curlH] d.xdr 

-2 J J (ar )~ 4 d a { [curl, u] ■ VH}(d Q curlH) dxdr 

+ 2 J J (ar)- 4 <9"curl((H ■ V)u - Hdivu) (S a curlH) dxdr. 
For | a | < 2, it is easy to see that 

f (ar)- 4 |a Q curlH| 2 dx < C* + tC{M). (2.33) 

\ a \<2 J 

Next, we treat the case \a\ = 3. By a straightforward calculation we arrive at 
((ar)- 4 d a curl H, <9 Q curl H) 

<C +tC(M) + 2 J J (ar)- 4 5 Q curl((H- V)u - Hdivu) (<9 Q curl H) da-dr 

=:C + tC(M) + 2 / I a (r)dT. 

Jo 

For the term I Q (r), using the fact that divH = 0, we have 

X] 7 «( r ) = X] / M- 4 d Q curl((H ■ v ) u - Hdivu)d Q curlHda; 

j Q |=3 M=3 

<-J2 J {(H" 4 5 Q (-(u-V)H + (H-V)u- Hdivu) 



|a|=3" 

x d a curl curl H} da; + C(M) 



< f{ar)- i d a ((H-\7)u-Hdivu)d a AUdx + C(M) 

\ a \=3 J 

< - J2 ! M~ 4 (H • V)<9 Q H9 Q (Au) dx 

- J (ar)" 4 (VA 2 HH)VAdivuda; + C(M), 
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where we have used the integration by parts several times in the last step. Thus 
the lemma is proved, taking into account the positivity of a(S, eq) and r(S, eq). □ 

Finally, puting Lemmas (|2.3[) - (|2.7[) together, and using the induction argument 
as in [25], we get the following estimate which finishes the proof of Theorem 1 2. II 

Lemma 2.8. There exist a constant Co > and an increasing function C(-) from 
[0, oo ) to [0, oo), such that for all e G (0, eo] and t G [0, T], 

||(5, q, u, H)(t)||2 < Co + (t + e)C(M(T)). (2.34) 

3. Incompressible limit 

In this section, we prove the convergence part of Theorem 11.11 by modifying 
the method developed by Metivier and Schochet [55], see also some extensions 
in UBS]. 

Proof of the convergence part of Theorem \l.l\ The uniform bound (|1.26[) implies, 
after extracting a subsequence, the following limit: 

(<? £ ,u £ ,H £ ) -» (g,v,H) weakly-* in L°°(0, T; H i {R d )). (3.1) 

The equations (fLlTf and (fl~T8]) imply that d t S e and <9 t H e G C([0, T},H 3 (R d )). 
Thus, after further extracting a subsequence, we obtain that, for all s' < 4, 

S^S strongly in C([0,T], H^(S. d )), (3.2) 

LT^H strongly in C([0, T], H^ c (R d )), (3.3) 

where the limit H G C([0, T], fff ' c (K d )) n L°°(0, T; H^ oc {R d )). Similarly, by (l2~22l 
and the uniform bound ()1.26|) . we have 

curl (r (S' e )u £ ) ^ curl (r (5)v) strongly in C([0, T], ff*'" 1 (M d )) (3.4) 

for all s' < 4, where r (S) = \im e ^ ro(S e ) := lim £ ^ r(<S' £ , 0). 

In order to obtain the limit system, we need to prove that the convergence in 
(|3Tj) holds in the strong topology of L 2 (0, T; #f ' c (R d )) for all s' < 4. To this end, 
we first show that q = and divv = 0. In fact, from (|2.14|) we get 

eE(S\ eg £ )9 t U £ + L(d x )U e = -eE(S c ,eq e )u e ■ VU £ + eJ £ . (3.5) 

Since 

E(S\eq t )-E(S%0)=O(e), 

we have 

eE{S\ 0)9 t U £ + L(9 X )U £ = eh £ , (3.6) 

where, by virtue of (fj~26l) . h £ is uniformly bounded in C{[0,T} 7 H 3 {R d )). Pass- 
ing to the weak limit in (|3.6|) . we obtain Vg = and divv = 0. Since q G 
I/°°(0,T;i? 4 (IR d )), we infer that q = 0. Now to complete the proof of Theorem 
11.11 we first prove the following proposition. 

Proposition 3.1. Suppose that the assumptions in Theorem \1.1\ hold, then q t 
converges strongly to in L 2 (0,T; Hi oc (M. d )) for all s' < 4, and divu £ converges 
strongly to in L 2 (Q,T; H^ 1 (R d )) for all s' < 4. 
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The proof of Proposition 13.11 is built on the the following dispersive estimates 
on the wave equations established by Metivier and Schochet [26] and reformulated 
in®. 

Lemma 3.2 ( [2j[26]). Let T > and uf be a bounded sequence in C([0,T], 
H 2 (R d )), such that 

e 2 d t (b £ d t w e ) - V ■ {c e Vw e ) = e £ , 

where e e converges to strongly in L 2 (0, T; L 2 (W d )). Assume in addition that, for 
some s > d/2+1, the coefficients (b e , c e ) are uniformly bounded in C([0, T], H 4 (M. d )) 
and converges in C([0, T], H^ oc (R d )) to a limit (b,c) satisfying the decay estimates 

\b(x,t)-b\ < Colx]- 1 - 6 , \V x b(x,t)\ < C \x\' 2 - S , 

\c(x,t) -c\ < Col.Tr 1 - 5 , \V x c(x,t)\ < C \x\- 2 - s , 

for some given positive constants b, c, Co and 5. Then the sequence w e converges 
toO in L 2 {Q,T;L 2 oc {R d )). 

Proof of Provosition \3.1[ We first show that q e converges strongly to in L 2 (0, T; 
Hi oc (R d )) for all s' < 4. An application of the operator e 2 dt to (|1.15|) gives 

e 2 d t (a(S\ eqt)d t q e ) + ed t divu e = ~e 2 d t {a{S\ e 9 e )(u e • V)q e }. (3.7) 

Dividing (|1.16|) by r e (S' e ,eg e ) and then taking the operator div to the resulting 
equations, one has 

<9 t divu e + idivf — — -Vq' 



= -div((u e ■ V)u e ) + div(-^— - -(curlH 6 ) x H e ). (3.8) 
Subtracting (|3.8[) from ()3.7|) . we get 

e 2 d t (a(ffVgW)-div( 1 w) = F(S e , g £ , u £ , H £ ), (3.9) 



where 

F(5 e ,q e ,u e ,H e ) =ediv( 1 6 . (curlH £ ) x H e 
\r[S e ,eq e ) 

- ediv((u e ■ V)u e ) - e 2 d t {a{S\ e<f)(u £ ■ V)g 6 }. 

In view of the uniform boundedness of (S e , q t 1 u e , H e ), the smoothness and pos- 
itivity assumptions on a(S e , eq e ) and r(S e , eg e ), and the convergence of S e , we find 
that 

F(S' e , g e ,u e ,H £ ) ->0 strongly in L 2 (0,T; L 2 (R d )), 



and the coefficients in (|3.9p satisfy the requirements in Lemma 13.21 Therefore, by 
virtue of Lemma [ 



q e ^0 strongly in L 2 (0, T; L 2 oc (R d )). 

On the other hand, the uniform boundedness of q e in C([0, T], if 4 (IR d )) and an 
interpolation argument yield that 

q^Q strongly in L 2 (Q,T; H^ c (R d )) for all s' < 4. 

Similarly, we can obtain the convergence of divw e . □ 
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We continue our proof of Theorem 11.11 From Proposition 13.11 we know that 

divu e ^divv in L 2 (0,T;H£- 1 (R d )). 

Hence, from (|3.4[) it follows that 

u e -> v in L 2 (0, T; H{ m (R d )) for all s' < 4. 

By (|3.2p . (|3.3p and Proposition ^. 2[ we obtain 

r e (S e ,eq e ) ->• r (5) in L°°(0, T; L°°(E d )); 
Vu £ ^ Vv in L 2 (0,T; i? 1 *'~ 1 (]R d )); 
VH £ ^ VH in L 2 (0,T;i7 1 s o ' ( T 1 (R d )). 

Passing to the limit in the equations for S e and H e , we see that the limits S and 
H satisfy 

d t S + (v • V)S = 0, <9 t H + (v • V)H — (H • V)v = 
in the sense of distributions. Since r(S e , eq e ) — r (S e ) = 0(e), we have 
(r(S\ eg") - r (S'))(d t u e + (u e • V)u e ) ^ 0, 

whence, 

r(S e ,eg £ )(d t u e + (u e • V)u e ) = (r(S e , e 9 £ ) - r (S £ ))(d t u e + (u e • V)u £ ) 

+ cMro(S e K) + (u e -V)MS e K) 
->m(S)(4v+(v-V)v) 

in the sense of distributions. 

Applying the operator curl to the momentum equations (|1.16[) and then taking 
to the limit, we conclude that 

curl (r (S)(d t v + v • Vv) - (curlfi) x H /jAv) = 0. 

Therefore, recalling curlV = 0, we see that the limit (S, v, H) satisfies 

r(S\0)(<9 t v+ (v • V)v) - (curlH) x H - ^Av + Vtt = 0, (3.10) 

<9 t H + (v - V)H - (H • V)v = 0, (3.11) 

d t S + (v • V)5 = 0, (3.12) 

divv = 0, divH = (3.13) 

for some function tt. 

If we employ the same arguments as in the proof of |26[ Theorem 1.5], we find 
that (S, v, H) satisfies the initial conditions (jl.28j) . Moreover, the standard iterative 
method shows that the system (|3.10[i — (|3.13|) with initial data f| 1 . 28[) has a unique 
solution (S*,v*,H*) £ C([0, T], iJ 4 (R d )). Thus, the uniqueness of solutions to the 
limit system f|3.10[) — (|3.13[) implies that the above convergence results hold for the 
full sequence (S e , q e , u e , H e ). Thus, the proof is completed. □ 
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